Antiferromagnetic exchange and spin-fluctuation pairing 
in cuprate superconductors 
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A microscopic theory of superconductivity is formulated within an effective p-d Hubbard model 
for a Cu02 plane. By applying the Mori-type projection technique, the Dyson equation is derived 
for the Green functions in terms of Hubbard operators. The antiferromagnetic exchange caused by 
interband hopping results in pairing of all carries in the conduction subband and high T c proportional 
to the Fermi energy. Kinematic interaction in intraband hopping is responsible for the conventional 
spin-fluctuation pairing. Numerical solution of the gap equation proves the d-wave gap symmetry 
and defines T c doping dependence. Oxygen isotope shift and pressure dependence of T c are also 
discussed. 



I. INTRODUCTION 

A unique property of cuprates is the huge antiferro- 
magnetic (AFM) superexchange interaction J ~ 1500 K 
which causes a long-range AFM order in the undoped 
regime, while in the metallic state results in strong AFM 
dynamical spin fluctuations. These are responsible for 
anomalous normal state properties of cuprates and for su- 
perconducting pairing as originally proposed by P.W. An- 
derson 1 . In a number of studies of the reduced one- 
band t-J model (see, e.g. 2 ' 3 and references therein) it 
was shown that the AFM exchange interaction results in 
ci-wave pairing with high T c . However, to asses the ori- 
gin of the AFM pairing mechanism one has to consider 
the original two-band p-d model for Cu02 plane without 
reducing it to an effective t-J model with instantaneous 
exchange interaction in one subband by projecting out 
the intersubband hopping. In fact, this transformation 
is used in derivation of the BCS model with an instan- 
taneous electron interaction by projecting out electron- 
phonon coupling in the original electron-ion model. 
' In this paper a microscopic theory of superconductiv- 
ity within the effective p-d Hubbard model is presented. 
By applying a projection technique for the matrix Green 
function in terms of Hubbard operators the Dyson equa- 
tion is derived 4 . It is found that the mean-field solution 
results in d-wave superconducting pairing mediated by 
the exchange interaction due to interband hopping, simi- 
lar to the t-J model. The self-energy caused by the kine- 
matic interaction in the intraband hopping is calculated 
in the non-crossing approximation. It defines the conven- 
tional spin-fluctuation pairing. Numerical solution of the 
gap equation gives the doping dependence of supercon- 
ducting T c and wave- vector dependence of the gap which 
are in agreement with experiments 4 . T c increase under 
pressure and oxygen isotope effect in cuprates are also 
discussed 3 . 



II. DYSON EQUATION 

Starting from the original two-band p-d model for 
Cu02 plane we can reduce it to an effective two-band 
Hubbard model by applying the cell-cluster perturbation 
theory (see, e.g. 5 ): 
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where X™ m = \in)(im\ are Hubbard operators for four 
states n,m = |0), \a), |2) = | fl), o = ±1/2, a = -a. 
Ei = ed — /j, and E% = 2E\ + A are energy lev- 
els of the lower Hubbard subband (one-hole Cu d-like 
states) and the upper Hubbard subband (two-hole p-d 
singlet states), respectively, /j, is the chemical potential, 
A = e p — €d ~ 3 eV is the charge transfer energy and 
t ~ A/2 is the p-d hybridization parameter. The hop- 
ping integrals | tf |< A, e.g. tf l±a = t nn ~ 0.14t - 
0.2 eV, which means that the Hubbard model (1) corre- 
sponds to a strong correlation limit. We treat it within 
the Hubbard operator technique to preserve rigorously 
a constraint of no double occupancy of any quantum 
state \in) which follows from the completeness relation: 
Xf° + Xf CT + Xf 5 + Xf = 1 . 

To discuss the superconducting pairing we introduce 
the four-component Nambu operators Xi„ and JtJ and 
define the 4x4 matrix Green function (GF) 6 : 
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Gij a and Fijcr are 



where x\ a = [Xf* Xf° Xf Xf) 
the two-subband normal and anomalous matrix compo- 
nents, respectively. By applying the projection technique 
for the equation of motion method for GF (2) we derive 
the Dyson equation in the (q, w)-representation : 
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where x = ({XioXj }) . The zero-order GF within the 
generalized mean field approximation (MFA) is defined 
by the frequency matrix which in the site representation 
reads 

E ija = AijaX' 1 , Aja = ({{X ia ,H],Xl}). (4) 

The self-energy operator in the Dyson equation (3) in the 
projection technique method is defined by a proper part 
(having no single zero-order GF) of the many-particle GF 
in the form 

E CT ( q ^) = rW r > \zn^ r 1 . (5) 

Here the irreducible Z-operator is given by the equa- 
tion: Z& r ^ = [Xi a ,H] — Y^i Euo-Xie which follows from 
the orthogonality condition: ({Z { J r) , x] a }) = 0. The 
equations (3)-(5) provide an exact representation for the 
GF (2). Its calculation, however, requires the use of some 
approximation for the self-energy matrix (5) which de- 
scribes finite lifetime effects (inelastic scattering of elec- 
trons on spin and charge fluctuations). 



III. MEAN-FIELD SOLUTION 

In the MFA the electronic spectrum and superconduct- 
ing pairing are described by the zero-order GF in Eq. (3) 
where the frequency matrix (4) reads: 
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Here u)ij a and Aij a are 2x2 matrices for the normal and 
anomalous components, respectively. The normal com- 
ponent defines quasiparticle spectra of the model in the 
normal state which have been studied in detail in 5 . The 
anomalous component defines the gap functions, e.g., for 
the singlet subband we have: 
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where Ni = J^a X™ + 2X 22 is the number operator. 
In terms of the Fermi operator: Ci a = X®' 7 + 2aX[ 2 
we can write the anomalous average as (cj^Cjf iVj) = 

(X° l Xj 2 Nj) = (XfNj) since other products of Hub- 
bard operators vanish according to multiplication rules 
for them: X^ 1 X^ = (5 7 ,a X" 13 . Therefore, the anoma- 
lous correlation functions describe the pairing at one lat- 
tice site but in different subbands. 

The same anomalous correlation functions were ob- 
tained in MFA for the original Hubbard model in 



Refs. 7-9 . To calculate them in 7 ' 9 the Roth procedure was 
applied based on decoupling of operators on the same 
lattice site in the time-dependent correlation function: 
{cii(t)\c^(t')Nj(t')) . However, the decoupling of Hub- 
bard operators on one lattice site is not unique and there- 
fore unreliable. The use of kinematic restrictions for the 
Hubbard operators 8 instead of decoupling also did not 
produce a unique solution. 

We calculated the correlation function (X® 2 Nj) with- 
out any decoupling from the equation of motion of the 
corresponding commutator GF L i3 \t — t') = ((Xf 2 (t) \ 
Nj(t'))) . By neglecting exponentially small terms of the 
order of exp(— A/T) we obtain for the correlation func- 
tion in the singlet subband (i ^ j) 4 : 



(X» 2 N j ) = -(l/A) 



2atli(X^X^N 3 ) 
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The approximate value is derived in the two-site approx- 
imation usually applied for the t-J model: X^ =J Nj = 
2X° 2 . Therefore, the gap equation (7) in the case of hole 



doping reads 



A 2 l = -2at\ 2 {X a t 2 N 3 ) = JijiXfX? 2 ), 
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where = 4(tj 2 ) 2 /A. It is the conventional exchange 
interaction pairing in the t-J model. For electron dop- 
ing analogous calculations for the anomalous correla- 
tion function ((2 — Nj)Xf 2 ) gives for the gap function 
A\j a = J tj {XfX«°). 

We conclude that the anomalous contributions to the 
zero-order GF, Eq. (3), are described by conventional 
pairs of quasi-particles in one subband and there are no 
new composite operator excitations, the "cexons", pro- 
posed in 9 . 



IV. SELF-ENERGY 

The self-energy matrix (5) can be written in the form 

(10) 
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where the 2x2 matrices M and $ denote the normal 
and anomalous contributions to the self-energy, respec- 
tively. We calculate them in the non-crossing, or the 
self-consistent Born approximation (SCBA). In SCBA, 
Fermi-like and Bose-likc excitations in the many-particle 
GF in (10) are assumed to propagate independently 
which is given by the decoupling of the corresponding 
time-dependent correlation functions as follows: 

{By (t)X 1 (t)B 2 , (t')X 2 (t')) ,2*2') 

^{ Xl {t)X 2 {t')){B v {t)B 2 ,{t')). (11) 
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Using the spectral representation for the correlation func- 
tions we get a closed system of equations for the GF (2) 
and the self-energy (10). For instance, the anomalous 
part of the self-energy for the singlet subband reads: 
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Here the spin-fluctuation pairing is defined by imaginary 
part of the dynamical spin susceptibility x'si'i^) = 
— (1/tt) Im((S 9 I S_«j)) w+i( 5 which comes from the cor- 
relation functions (Bi>(t)B 2 i(t')) . 

For the hole doped case, /x ~ A, the poles of the GF 
F% 2 (k, uji) are close to the Fermi surface (FS), while the 
poles of the GF F ( J 1 (k, cji) are situated much below the 
FS at u>i ~ A 3> W where W is an effective bandwidth. 
Therefore, we can neglect the second term and use the 
weak coupling approximation (WCA) for the first term 
in the self-energy (12). By taking into account the ex- 
change interaction contribution (9), we get the following 
equation for the singlet gap ^(q): 



*(q) = ^E^/( k ' < l) 
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interaction 14/ (k, q) = J(k — q) — 
q, u = 0). The quasiparticle energy in the 



where the 
*22(k)x s (k 

singlet subband is given by E 2 (k) = [fi 2 (k) + $ 2 (k)] x / 2 
where f2 2 (k) is the energy in the normal state 5 . For an 
electron doped system (fi ~ 0) the WCA equation for 
the gap <& n (q) is similar to Eq. (13). 



V. RESULTS AND DISCUSSION 

Let us compare the exchange and spin-fluctuation pair- 
ing in the gap equation (13). The exchange interac- 
tion for nearest neighbors J(q) = 4J~/(q), j(q) — 
(1/2) (cos q x + cosq y ) . The spin-fluctuation contribu- 
tion is defined by hopping interaction f 22 (k) ~ t„„7(k) 
with t nn ~ 0.2 eV and the static susceptibility which we 
model by the function: x«(q) = XoM(l + £ 2 [1 + 7(q)]), 
where £ is the AFM correlation length, u> s ~ J, and xo 
is defined from the normalization condition: (SiSi) = 
E q (S q S q ) - (3/4)(l -6),S=\l-n\. 

By adopting a model for the d-wave gap function 
$(q) = (f>(cosq x — cosq y ) = <f>r](q) we get from (13) the 
following equation for T c : 
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where X s ~ t 2 in juj s . For the exchange interaction medi- 
ated by the interband hopping with a large energy trans- 
fer A » W the retardation effects are negligible that 



results in pairing of all charge carriers in the conduction 
subband: 



VlMW - /u)exp(-l/A ex ), 



(15) 



where X ex ~ J N(p) , N(^i) is the density of electronic 
states. For spin-fluctuation pairing we obtain a conven- 
tional BCS-likc formula for T r : 



T c ~ w s exp(-l/A s/ ) . 



(16) 



where A s / ~ X s N(p). By taking into account both 
contributions, a result similar to (16) is obtained, 

with an effective coupling constant: = + 

A e2; /[1 — X ex ln(^/w s )] . By taking for estimate \i = 
W/2 ~ 0.35 eV, uj s ~ J ~ 0.13 eV and A s / ~ X ex = 0.2 

we get x[ e / f) ~ 0.2 + 0.25 = 0.45 and T c ~ 160 K, 
while the spin-fluctuation pairing alone gives T° ~ 
Lo s exp(— 1/A S /) ~ 10 K. 

A direct numerical solution of the gap equation (13) 
in 4 proved the d-wave symmetry of the gap and gave 
the maximum T c ~ 0.12t rm for /1 ~ W/2 at an optimal 
doping 8 pt — 0.12. The spin- fluctuation interaction pro- 
duces much lower T c since it couples the holes in a narrow 
energy shell, uj s <C /1 , near the FS. Moreover, it vanishes 
along the lines \k x \ + \k y \ — it where 7(k) = 0. 

While in clectron-phonon superconductors T c de- 
creases under pressure, in cuprates it increases. In par- 
ticular, in mercury superconductors dT c /da ~ —1.35 x 
10 3 K/A 10 (a is the in-plane lattice constant) and 
d\nT c /d\na ~ —50. From Eq. (15) we get an estimate: 



d\nT r d\nT r d\nj 



d\na 



d\n,J din < 



14 
T 



-45, 



(17) 



for X ex ~ 0.3, which is quite close to the experimentally 
observed one. Here we take into account that J(a) oc t^ d 

and t p d(a) oc l/(a) 7 / 2 . 

To explain a small oxygen isotope effect in cuprates, 
a c = -dlnT c /dlnM < 0.1, we can use Eq. (15). By 
taking into account the experimentally observed isotope 
shift for the Neel temperature in La 2 Cu04 11 : c*n = 
-d\uT N /d\nM ~ 0.05 with T N oc J we get 



a c = - 



rflnT c 
rflnM 



dlnT, dlnT N 



d\n,J d\nM 



0.16, 



(18) 



for A ~ 0.3 which is close to experiments. 

To conclude, the present investigation proves the ex- 
istence of a singlet c? a ,2_ y 2-wave superconducting pairing 
for holes or electrons in the p— d Hubbard model. It is me- 
diated by the exchange interaction and antiferromagnctic 
spin-fluctuation scattering induced by the kinematic in- 
teraction, characteristic to the Hubbard model. These 
mechanisms of superconducting pairing are absent in the 
fermionic models (for a discussion, see Anderson 12 ) and 
are specific for cuprates. 
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